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ABSTRACT 

We report the discovery of a medium-strength (~0.5 kG) magnetic field on the 
young, massive star r Sco (BO. 2 V), which becomes the third- hottest magnetic star 
known. Circularly polarized Zeeman signatures are clearly detected in observations col- 
lected mostly with the ESPaDOnS spectropolarimeter, recently installed on the 3.6-m 
Canada-France-Hawaii Telescope; temporal variability is also clearly established in 
the polarimetry, and can be unambiguously attributed to rotational modulation with 
a period close to 41 d. Archival UV spectra confirm that this modulation repeats over 
timescales of decades, and refine the rotation period to 41.033 ± 0.002 d. 

Despite the slow rotation rate of r Sco, we nonetheless succeed in reconstructing 
the large-scale structure of its magnetic topology. We find that the magnetic structure 
is unusually complex for a hot star, with significant power in spherical-harmonic modes 
of degree up to 5. The surface topology is dominated by a potential field, although a 
moderate toroidal component is probably present. We fail to detect intrinsic temporal 
variability of the magnetic structure over the 1.5-yr period of our spectropolarimetric 
observations (in agreement with the stable temporal variations of the UV spectra), 
and infer that any differential surface rotation must be very small. 

The topology of the extended magnetic field that we derive from the photospheric 
magnetic maps is also more complex than a global dipole, and features in particular 
a significantly warped torus of closed magnetic loops encircling the star (tilted at 
about 90° to the rotation axis), with additional, smaller, networks of closed field lines. 
This topology appears to be consistent with the exceptional X-ray properties of r Sco 
and also provides a natural explanation of the variability observed in wind-formed UV 
lines. Although we cannot completely rule out the possibility that the field is produced 
through dynamo processes of an exotic kind, we conclude that its magnetic field is 
most probably a fossil remnant from the star-formation stage. 
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1 INTRODUCTION 

Magnetic fields in hot, high-mass stars of spectral types O 
and early B may have a drastic influ ence on the physics of 
the stellar interiors fe.g..'S Druitl200 and atmospheres (e.g., 
Babel & Montmerle 1997). As a consequence, they can also 
significantly modify these stars' lon^-ter m evolution, and in 
parti c ular their ro tational history (e.g., iMaeder Mevnetl 
I2OOI l2004 120051) . However, quantifying these effects re- 
quires that we know the basic properties of such fields. 

From a theoretical point of view, several quite differ- 
ent mechanisms have been proposed as potential means 
of generating large-scale magnetic fields in very hot stars. 
For example, such fields may be fossil remnants of the 
star-formation stage, either as relics of the field that per- 
vaded the interstellar medium from which the star formed, 
or as leftovers from dynamo action in the convective 
Hayashi phase. This hypothesis was initially proposed for 
the lower- mass, chemic ally peculiar, magnetic Ap and Bp 
stars (e.g.. Moss 2001) and has recently been rediscussed 
for the particular case of high- mas s O an d early- B stars by 
iFerrario Wickramasingha (|2005L 120061) : in this scenario, 
very hot magnetic stars are expected to be the progenitors 
of magnetic neutron stars. 

A second possibility is that fields may be pro- 
duced by continuing dynamo action; this optio n, ini- 
tially suggested about two decades ago ("e.g.. IMos3|1982|) . 
has not been closely examined until quite recently, 
but is now attracting increasing interest from theoreti- 
cians. A dynamo action could opera t e in t he convec- 
tive core ((Charbonneau MacGrego3 1200 ll : iBrun et all 
120051 ). but the outer radiative envelope could also be 
involved in genuine dyn amo process es, in a subsurface 
shear layer (|Tout Pring le 1995; Lignieres et al. 1 996 ) 
or even throughout the whole envelope (ISpruitl 
MacDonald Mullan |_[2004l: iMullan MacDonaldl E 
Maeder M evnet 2OO.5I: lBraithwaitell2006^ ! 



Unfortunately, all available theoretical options still suf- 
fer significant problems. Fossil- field theories have yet to 
demonstrate that sufficient magnetic flux can survive the 
accelerated decay and expulsion associated with Hayashi 
convection^, while core-dynamo theories are still lacking de- 
tailed models linking the field produced in the core with that 
emerging at the surface ( Moss 200 1). Being based on a very 
new idea, radiative-zone dynamo theories are still in their 
infancy and need to establish their potential validity, by 
demonstrating that, in particular, stellar radiative zones are 
capable of sustaining differential rotation throughout their 
whole volume. 

In addition to these issues, existing theories face 
a number of problems when compared to observations. 
At some point, all theories involve the coexistence of 
large-scale toroidal fields in st ellar interiors and possibly 
even close to the surface (e.g. , iBraithwaite Spruitll200'3: 
iBraithwaite NordhlIldll2006^ . which o bservations of mag- 
netic Ap and Bp stars do not yet confirm (|MosJ200 J) . More- 
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^ This difficulty may not be as problematical as it first seemed, 
as discussed in Moss (2003), and should not concern stars with 
masses larger than 10 Mq. 



over, dynamo theories predict magnetic topologies that are 
expected to vary on short timescales and to depend strongly 
on stellar rotation rates, which again is not observed in 
intermediate-mass magnetic stars (Moss 2001). Last but not 
least, dynamos theories should in principle apply to most hot 
stars (as they do for most cool stars), making it hard to un- 
derstand why only a small fraction of them (about 10% in 
the case of Ap/Bp stars) apparently hosts magnetic fields. 



Indirect evidence for the presence of magnetic fields in 
high- mass stars is regularly reported in the literature, these 
being postulated as a potential explanation for many oth- 
erwise enigmatic phenomena, including unantic ipated X-ray 
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line profiles and high X-ray temperatures (e.g., 
'1997; "Robinso n etllll2002l: ICohen et al.1 1200^1 : 

2004). However, with di rect detections of m agnetic fields in 
only two O stars to date (iDonati et al. ''2002. 2006), together 
with less than a ha ndful of early-B sta r s with masses large r 
than 10 Mo (e.g., IPonati et al.ll200ll : iNeiner et al.ll200,^) . 
very little is known reliably about magnetic strengths and 
topologies from a statistical point of view. One reason for 
this is that absorption lines of these stars are not only rel- 
atively few in number in the optical, but are also generally 
rather broad (due to rotation, or to some o ther type of as 
yet u nknown macroscopic mechanism; e.g., 'Howarth e t al.l 
Il997r l. decreasing dramatically the size of the Zeeman sig- 
natures that their putative fields can induce. 



In these respects, r Sco (HR 6165, HD 149438, 
HIP 81266) appears as a very promising target for magnetic- 
field studies; thanks to its brightness and unusually nar- 
row absorption lines (among the sharpest known for stars 
more massive than 10 M©), r Sco is an obvious candidate 
for accurate spectropolarimetric experiments. Moreover, its 
strong, hard X-ra y emission (with logLx/i^Boi — —6.5, 
ICohen et al1ll997l) ' poses a severe challenge to the standard 
picture of 0-star w ind-shock models, leading some authors 
(e.g., Cohen et aLl boO^") to speculate that it displays the 
presence of a mag netic ally confined wind, a s prop osed by 
Iptahl et al.| ^22^ and|Babel <^ Montmerlej (Il 99/t. and as 
detected on the very you ng O star 0^ Ori C (|Donati et al.l 
l2002l :lrTagne et al.' 2005b ?V We therefore included r Sco in 
a list of candidates for observation with ESPaDOnS, a high- 
efficiency spectropolarimeter installed in 2004/5 on the 3.6- 
m Canada-France-Hawaii Telescope (CFHT; Donati et al., 
2006, in preparation). 



In this paper, we first present our new spectro- 
polarimetric observations and Zeeman detections obtained 
for T Sco, and briefly summarize the UV spectra (Sec. 0). 
We then carry out a basic analysis of these data to establish 
the rotation period (Sec. |3, and review the fundamental 
stellar parameters in the light of our result (Sec^J. We per- 
form detailed modelling of the magnetic topology of r Sco, 
by direct fitting to the observed Zeeman signatures (SecEJ. 
Finally, we discuss the implications of our results for models 
of the X-ray emitting magnetosphere of r Sco, as well for as 
theories of large-scale magnetic- field generation in high- mass 
stars (Sees. island 0. 
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Table 1. Journal of spectropolarimetric observations. Columns 1-6 list the UT date Sz time, heliocentric Julian date (all at mid-exposure), 
observing site, exposure time, and peak signal to noise ratio (per 2.6 kms"-*^ velocity bin) of each observation. Columns 7 and 8 list the 
rms noise level (per 1.8 kms"-*^ velocity bin, relative to the unpolarized continuum level Ic) in the circular-polarization profile produced 
by Least-Squares Deconvolution (see Sec.|2| and the estimated longitudinal field (with corresponding Icr error bars). The rotational 
cycle E from the ephemeris of ean. ITlis given in column 9. 
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2 OBSERVATIONS 

2.1 Optical spectropolarimetry 

Spectropolarimetric observations of r Sco were collected 
with ESPaDOnS from 2004 Sep. to 2006 Feb., during the 
first engineering runs (in 2004, when the field was first 
detected) and subsequently for scheduled ESPaDOnS pro- 
grammes. The ESPaDOnS spectra span the entire optical 
domain (from 370 to 1,000 nm) at a resolving power of about 
65,000. In total, 32 circular-polarization sequences were col- 
lected, most of them consisting of 4 individual subexposures 
taken in different polarimeter configurations (in two cases, 
the sequence was interrupted after the second exposure due 
to a technical problem with the in strument con t rol). The 
extraction procedure is described in lDonati et al.l (|l997l see 



Donati et al., 2006, in prep., for further details) and was 
carried out with Libre ESpRIT, a fully automatic reduction 
package/pipeline installed at CFHT for optimal extraction 
of ESPaDOnS spectra. The peak signal-to-noise ratios per 
2.6 kms~^ velocity bin range from 490 to 1710, depending 
mostly on weather conditions and exposure time (see Ta- 
blelD^. 

Three additional circular-polarization spectra were ob- 
tained at the 3.9-m Anglo- Australian Telescope (AAT) in 
late Sep. 2004, using a visitor-instrument Cassegrain po- 

^ Note that the instrument suffered a 1.3-mag loss in throughput 
compared to the optimal performance from early March to late 
June 2005, due to severe damage to the external jacket of optical 
fibres linking the polarimeter with the spectrograph (now fixed). 
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LSD profiles, r Sco, ESPaDOnS, 2005 Jun. 24 




-50 



Velocity (km/s) 

Figure 1. LSD unpolarized and circularly polarized profiles of 
r Sco (bottom, top curves respectively), as observed on 2005 
June 24. The mean polarization profile is expanded by a factor of 
100 and shifted upwards by 1.10 for display purposes. 



larimeter, fibre linked to the UCL echelle spectrograph 
(UCLES), in a setup very similar to that described by 
iDonati et all ((20 03 "l. The A AT spectra, processed with the 
same reduction package as used for the CFHT data, cover 
430-670 nm at a resolving power of 65,000. The peak signal- 
to- noise ratios per 2.6 kms~^ velocity bin range from 670 to 
890, depending on the weather (see Table CJ. 

Both instruments use Fresnel rhombs (rather than crys- 
talline plates) as retar ders, with the result that spectropo - 
larimetric ripples (e.g.. lAitken HoughlboOll : ISemeJl2003h 
are decreased down to a level below detectability. An ex- 
ample of this improvement compared to previous otherwise 
sim ilar instruments (su ch as the MuSiC oS spectropo l arime- 
ter: jPonati et al.lll99ftl is illustrated bv lWade et all (|200i) 
and' W ade et all (|2006 h in the particular case of hot-star ob- 
servations. 

In order to gain a multiplexing advantage by combin- 
ing results from different lines in the spe ctra, Least- Squares 
Deconvolution (LSD; Don ati et al.lll99'^ was applied to all 
observations. The line list required for L SD was c omputed 
from an Atlas9 LTE model atmosphere (|Kuruc7l [l993^ at 
Teff = 30,000 K, \ogg = 4.5, roughly matching r Sco's pa- 
rameters. We utilized only moderately strong lines (those 
with synthetic profiles having line-to-continuum core depres- 
sions larger than 10% prior to all non-thermal broadening 
mechanisms, but omitting the very strongest, broadest fea- 
tures, such as Balmer and He lines, whose Zeeman signature 
is strongly smeared out compared to those of narrow lines) 
- some 500 spectral features altogether, with about half cor- 
responding to oxygen lines and most others coming from N, 
Si, C and Fe. The average noise levels in the resulting LSD 
signatures range 0.6-2.1x10"^ for ESPaDOnS spectra and 
1.0-1.4x10"^ for the AAT data (per 1.8-kms~^ velocity bin, 
relative to the unpolarized continuum level Ic] see Table CJ. 
Significant Zeeman signatures are clearly detected in all the 
spectra, with a full amplitudes of about 0.2%, demonstrat- 
ing that a magnetic field is securely detected for r Sco; an 
illustrative example is shown in Fig. Q 



The Zeeman signatures detected are variable with time^ 
and correspond to projected longitudinal fields (computed 
from the f irst-order moment of the Stokes V LSD profile; 
iDonati et al . 1997) ranging from -50 to +90 G (see Ta- 
ble^ -^Our results are compatible with the earlier findings of 
uLandstreetl (119821 ) that longitudinal fields of normal upper- 
main- sequen£e_staTsJiigeneral, and of r Sco in particular 
(for which Landstr eetl[l982 report two estimates), are less 
than 100 G on average. 

2.2 UV spectroscopy 

We augmented the new spectropolarimetric data set with 
archival UV spectra obtained with the International Ultra- 
violet Explorer satellite (lUE; Boeeess et al. 1978). The lUE 
archive contains 107 high-resolution spectra of r Sco, cov- 
ering the wavelength range 115-195 nm, sampled at 0.01- 
nm intervals with a resolution of about 10^ and S/N ~ 
20. We examined spectra obtained from the MAST lUE 
archive at the Space Telescope Science Institute, reduced us- 
ing NEWSIPS pipeUne software ((Nichols Linskvlll996l the 
'mxhi' product). lUE spectra were each acquired through 
one of two focal-plane apertures, designated 'small' and 
'large'; we found indications of data-processing errors in a 
number of small-aperture spectra, which also generally have 
poorer signal and S/N than the large- aperture data (typi- 
cally, the small aperture transmitted only ^half the incident 
radiation). We therefore present results based on only the 
72 large- aperture spectra in the archive; incorporating the 
small-aperture data introduces no material changes to the 
conclusions. 



3 THE ROTATION PERIOD OF r SCO 
3.1 Zeeman spectroscopy 

Time variability of the projected longitudinal field is a po- 
tentially powerful tool for estimating stellar rotation peri- 
ods, which often cannot be measured by other means. Specif- 
ically, temporal fluctuations of the line-of-sight component 
of the field in hot stars may be attributed to a magnetic 
topology that is not axisymmetric about the rotation axis, 
thus showing different configurations to the observer as the 
star rotates. Searching for timescales on which the longitu- 
dinal field repeats identically from one cycle to the next has 
been very successf ul in estimating; rotation periods of Ap 
and Bp stars (e.g.,|Borra <^ Landstree^ll^Sd) • 

Of course, the longitudinal-field values contain far less 
information about the field topology than the Zeeman sig- 
natures themselves. The Stokes V LSD profile shown in 
Fig.Qprovides a clear demonstration of this point; while the 
Zeeman signature strongly indicates a field detection (at a 
level of about 27cr), the corresponding longitudinal field is 
1.3 =b 3.1 G and is therefore inconclusive by itself. Nonethe- 
less, longitudinal-field values are convenient summary statis- 
tics in several respects, and in particular can provide useful 

^ The LSD / profiles are also variable in strength. The level of 
this variability is, however, small (about 0.3% rms of the contin- 
uum level, equivalent to ~3% of the central line depth), affecting 
mainly the line width. We return to this point in Sec. ^ 
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rough estimates of the average magnetic flux over the stehar 
surface. 

As an initial step in modelling the longitudinal-field 
variations of r Sco, we compared a double sine-wave fit to 
the data, with the two periods held in the ratio 2:1 (moti- 
vated by what we expect from a simple, rotationally modu- 
lated, linear combination of dipole plus quadrupole fields). 
For a range of assumed values for the longer (i.e., rotational) 
period, the amplitudes and phases of both waves were opti- 
mized to obtain the best fit to the longitudinal-field measure- 
ments, using a standard least-squares minimization process, 
with the reduced-x^ value, x^i evaluated as the statistic of 
merit for the fit quality. Fig. 12 (upper panel) shows the re- 
sults of this exercise, and indicates a strong minimum in 
at Prot = 41.08 ± 0.07 d (la uncertainty). 

However, the double sine- wave model evidently provides 
no more than a rough fit to the data, with a minimum 
as large as '-^9. This, in turn, indicates that the field 
must be significantly more complex than a tilted dipole plus 
quadrupole combination, and must contain higher-order, 
multipole components.^ Clearly, more-detailed modelling of 
the Zeeman signatures is necessary in order to obtain an 
adequate description of the magnetic topology of r Sco; we 
describe such modelling in Sec. El 

The mis-match between the simple model and the data 
means that the 41-d signal is no more than a provisional 
estimate of the rotation period. More accurate and precise 
values may be obtained both from modelling the Zeeman 
signatures directly (Sec.|5|, and from a time-series analysis 
of the lUE data (Secg^. 



3.2 UV line-profile variability 

Direct measurement of magnetic fields in hot stars (e.g.; 
i Donati et all l200lL 120021: iNeiner et all 1200.4 IPonati et all 
I2OO6I) has often been presaged by the detection of strictly 
periodic UV line-profile variability, with the inference of 
a magnetic field by analogy with He peculiar stars (e.g.; 
IStahl et al.lll996l : iNeiner et all 1200^1: IWalborn e t al. 20041). 
A search for variability in the lUE spectra of r Sco 
(known to exhibit abn ormally strong UV P-Cygni lines; 
IWalborn Paneklll984ri is therefore a natural follow-up to 
our detection of rotationally modulated Zeeman signatures 
in the spectropolarimetric data. 

We are fortunate that r Sco was adopted as a photo- 
metric calibration star by the lUE Project, and that as a 
consequence the temporal sampling of the UV spectra is 
rather well suited to investigating the rotational timescale 
identified in the Zeeman spectroscopy: the large- aperture 
spectra span 16.6 yr (1979 Feb - 1995 Sept, median date 
1988 Aug) at a median sampling rate of 50.7 d (range 0.7 hr 
-417 d). 

We analyzed the rectified lUE spectra using the 
implementation of the clean algorithm described by 



The rotational modulation of the longitudinal field is very rem- 
iniscent of what is observed on the helium-strong star HD 37776 
jThompson Landstreetll98f^. apart from the fact that the field 
of HD 37776 is stronger than that of r Sco by more than an order 
of magnitude. 



Periodogrann fronn longitudinal field data 
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Figure 2. Periodogram resulting from a double sine- wave fit to 
the longitudinal-field data of r Sco. Top panel: xi cis a function 
of the period of the main sine wave. A clear minimum is obtained 
for a period of about 41 d. Middle panel: Temporal fluctuations of 
the longitudinal field of r Sco (full dots, with la error bars) along 
with the model fit (full line) for the adopted period of 41.033 d, as 
a function of heliocentric Julian date (HJD). Bottom panel: As for 
the middle panel, but as a function of rotation phase, computed 
using the ephemeris of eqn. Q 
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Figure 3. Time-series analysis of lUE spectra of r Sco. The bottom panel shows the mean spectrum for (left to right) the N v, Si IV, and 
C IV resonance doublets; in each case, vertical tickmarks indicate the rest wavelengths, while the horizontal bar extends from 1000 kms"-*^ 
bluewards of the short-wavelength component to +500 kms~^ redwards of the long- wavelength component. The next panel up shows 
part of the 2D CLEANed periodogram, showing evidence for a periodic signal at u = 0.049 d"-*^ in the wind-formed lines (see Fig. 
while the next two panels show the fourier power and phase, respectively, at that frequency. The top panel shows the residuals about 
the mean of the rectified spectra phased according to the ephemeris of eqn. 
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Figure 4. Mean power spectrum for regions extending 
1000 kms"-*^ bluewards of the short-wavelength component to 
+500 kms"-*^ redwards of the long- wavelength component for the 
N V 124 nm. Si IV 140 nm, and C iv 155 nm resonance doublets 
(black line), and for equivalent regions 2 nm longwards (red line). 



iR.oberts et alJ [l98?); some results of this analysis are pre- 
sented in Fig. El There is a strong signal v — 0.049 
(P = 20.5 d) in the main features formed in the stellar wind 
(the N V 124 nm, Si IV 140 nm, and C IV 155 nm resonance 
lines, but also Si ill 120.6 nm). To examine this in greater 
detail, we constructed the mean power spectrum for wave- 
lengths from 1000 kms~^ bluewards of the short- wavelength 
component to 500 kms~^ redwards of the long- wavelength 
component for each of the three resonance doublets, together 
with a comparison spectrum from equivalent regions 2 nm 
longwards of each feature; results are shown in Fig. 2| 

In addition to the main signal, there is a weak secondary 
signal dX v — 0.024 d~^ (Fig. El? corresponding to a pe- 
riod of about 41 d.^ Since a period of 20.5 d is completely 
excluded by the Zeeman measurements (Sec. 13. If) , we in- 

^ There are also signals, in both the main and reference datasets, 
at P = 0.5 yr and 1 yr. These are presumed to be non- 
astrophysical (and are most probably the result of varying back- 
ground levels). 
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Figure 5. Mean residual intensity for regions from 1000 kms~^ 
bluewards of the short-wavelength component to +500 kms"-*^ 
redwards of the long- wavelength component for each of the 
N V 124 nm, Si IV 140 nm, and C IV 155 nm resonance doublets, 
as a function of rotational phase. Phases are computed according 
to the ephemeris of eqn. 



terpret our findings as indicating a slightly non-sinusoidal 
double wave with Prot = 41.03 d. To refine this period we 
average results from the 75 0.1- A wavelength bins in the 
— 1000/ + 500 kms~^ test ranges for which the peak of the 
power spectrum is safely identifiable with the f = 0.049 d~^ 
signal; half the mean peak frequency for those samples cor- 
responds to Prot = 41.033 ± 0.002 d (s.e.). 

As an alternative characterization of the data, fit- 
ting a gaussian to the peak of the mean power spectrum 
yields Prot = 41.034 d, and a full- width at half-maximum 
(FWHM) of 0.36 d. To obtain a second, reasonably con- 
servative, estimate of the uncertainty on the period, we may 
assume that the signal does not go out of phase by more than 
0.125Prot (i.e., 0.25 of the 20.5-d half-wave period) in each of 
the contributing wavelengths, giving P = 41. 0340 ±0.0040 d 
(s.e.). 

We emphasize that the data are fully consistent with 
a strictly periodic signal; the lUE data, spanning almost 
17 yr, show no evidence for any non-periodic component 
insofar as the width of the periodogram peak is entirely ac- 
counted for by the finite number of cycles in the dataset 
(iV = 148). Moreover, the same period is recovered, to high 
accuracy, in completely independent datasets separated by 
a quarter- century (viz., the lUE and Zeeman spectroscopy; 
more-detailed modelling of the Zeeman data, reported in 
Sec.lSl yields Prot = 41.02 ±0.03 d). Furthermore, the preci- 
sion of the adopted rotation period is high enough to ensure 
that the relative phasing between the archival lUE spectra 
and our new spectropolarimetric data is better than 1%. 

On the basis of these results we therefore adopt 



To = HJD 2, 453, 193.0 ± 41.033(±0.002)^ 



(1) 



as our rotational ephemeris, where phase zero is arbitrarily 
chosen for convenience as a date just prior to the acquisition 
of the first Zeeman data. The rotational modulation of the 
mean residual intensity over the N v. Si IV and C IV UV 
resonance doublets with this ephemeris is shown in Fig. |5| 
significant absorption episodes are observed to occur in all 
3 spectral features around phase 0.3 and 0.8. 



4 ROTATIONAL PROPERTIES 

From our analysis, r Sco is the second- slowest rotator so far 
known among high-mass stars (cp. HP 191612, whose r ota- 
tion period is suggested to be 538 d:'Don ati etaLll200(t) . In 
order to be able to refine our modelling of its magnetic topol- 
ogy, we first summarize what we know about its physical 
parameters, and examine how these tie in with the rotation 
period we derive. 

From deta iled spectroscopic studies by Kilian fT992^, 
iMokiem et ahh ^^OOS") and Simon -Diaz et al. (2006), we know 
that T Sco is a young star, with an age of a few Myr^, a mass 
of ~15 M0, a temperature of 31, 500 ± 500 K, and a radius 
of 5.2 ± 0.5 R©/ The radius and rotation period imply an 
equatorial rotation velocity, Ve^ of only 6.4 ± 0.6 kms~^, 
and thus a Une-of-sight projected rotation velocity, Ve sin(i), 
which is even less (where i is the inclination of the rota- 
tion axis to the line of sight). This Vesm.(i) lim it is s ignifi- 
cantly lower than the estimate given by Kiliaij (|l992l) from 
line-profile modelling (19 kms~ ^), but is simila r the value 
(of 5 kms~M report ed by both ISmith Kard ([l978 ^ and 
iMokiem e^_a]jj|2 005f) . and is cons istent with the upper limit 
given bv 'Simon- Di az et al.l (|2006l . 13 kms~^). 

We performed our own simple modelling of the LSD 
Stokes / profile, to get an estimate of how much (and what 
type of) broadening is required in addition to a pure ro- 
tational broadening in order to obtain a reasonable fit to 
the data. We find that gaussian broadening with an average 
FWHM of 13.5 kms~^ gives a nice match to the mean LSD 
profile (after allowance for a gaussian instrumental broad- 
ening, with a FWHM of 5 kms~^). This corresponds to a 
turbulent veloc ity of about 8 kms~\ a va lue similar to esti- 
mates given by iMokiem et al.l (|2p05") and Simon- D iaz et aTl 
(2006, 10.8 and 8.7 kms~^, respectively). This broadening is 
larger than the thermal broadening expected for most lines 
used in the LSD process (about 5 kms~^ in average, ranging 
from 3 to 6 kms~^, depending on the atom), and the similar 
widths found for lines from different atoms confirm that the 
main origin is non-thermal. 

As already mentioned in Sec. |21 we observe that the 
Stokes / LSD profiles of r Sco are slightly variable with time, 
and that this temporal variation can also be attributed as 
modulation with the 41 d rotation period derived at Sec. 
(see Fig. This modulation consists mainly of a small 
^10%), phase-locked variation of the line width, with the 
average photospheric profile being slightly narrower around 
phase 0.5. We suspect that this indicates that the atmo- 
spheric turbulence evinced by the line broadening is induced 
at the base of the wind, and is not constant over the surface 
of r Sco, reflecting, perhaps, the spatial distribution of field 
strengths and orientations; another possibility is that the 
density at the base of the wind is not constant over the sur- 
face of the star (again, as a potential result of the magnetic 



^ Alth o ugh T Sco is quoted to be younger than 1 Myr by 
|Kilian| ll99^, its m embership in the Sco OB2 association 
Jde Zeeuw et al.lll999^ suggests that r Sco, as the most massive 
star in the group, is roughly as old as the group itself - i.e., nearly 
5 Myr. 

^ Significantly l arger radii are cited in the literature, e.g. by 
iHowk et al.l l200n. but they are not consistent with the well- 
established Hipparcos parallax. 
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Figure 6. Dynamic LSD Stokes / spectrum of r Sco, demonstrat- 
ing that the average unpolarised hne is also modulated with the 
detected rotation period of 41 d (see Sec^J- A simple synthetic 
profile (featuring gaussian turbulence broadening with a FWHM 
of 13.5 kms"-*^ and VeS\n{i) = 6 kms"-*^, see text) was removed 
from all profiles to emphasize variability. A clear signal, centred 
at phase 0.5, crosses the line profile and indicates that the aver- 
age photospheric line of r Sco gets slightly narrower (by about 
10%) around phase 0.5. Blue and red respectively correspond to 
line absorption/emission features with an amplitude of zbO.4% of 
the continuum intensity (i.e. ±3% of the central line depth). 



field), causing the observer to see to different atmospheric 
depths (and thus different turbulent broadenings) at differ- 
ent rotational phases. We also find that a slightly better 
fit to the Stokes V signatures is obtained when assuming a 
local line profile width about 10% narrower than that de- 
rived from LSD Stokes / profiles; we speculate that Zeeman 
signatures are formed at atmospheric depths slightly differ- 
ent from the unpolarised spectral lines, i.e. in a layer where 
turbulent broadening is weaker. 

From the observed variability of the longitudinal field, 
which in particular shows a steep gradient as well as a sign 
switch in as little as 10% of the rotational cycle (between 
phase 0.85 and 0.95; see Fig.|21 bottom panel), we conclude 
that it is unlikely that r Sco is seen near-pole-on; modelling 
of the Zeeman signatures (Sec. reveals that significantly 
better fits to the data are obtained with Ve sin(i) = 6 km s~"^ 
than with values <5 kms~^. This implies that i is large, 
probably in the range 60°-90°. Since there is only a low 
probability that the star is seen exactly equator on, we as- 



sume i = 70° (and Veshi{i) = 6 kms~^) in the following 
analysis. Very similar results are obtained with values of i 
ranging from 60° to 90°. 

Finally, we note that the radial velocity of r Sco has 
remained remarkably stable during the whole length of our 
ru n, at — 0.6di0.1 kms~^. T his confirms earlier conclusions 
bv lStickland Llovdl ()l995f) that there is no reason to sup- 
pose that T Sco a member of a binary system. 



5 MODELLING THE SURFACE MAGNETIC 
TOPOLOGY OF r SCO 

5.1 Methodology 

To reconstruct the surface magnetic topology of r Sco from 
the set of observed Zeeman signatures, we use our m agnetic- 
imaging code ( Brown et al.iri991 ; Donati Brownlil997) in 
its most recent implementation (|Donatill200lh . While still 
based on maximum-entropy image reconstruction, this lat- 
est version reconstructs the field topology as a spherical- 
harmonic decomposition, rather than as a series of indepen- 
dent magnetic- image pixels as before. One obvious advan- 
tage of this method is that we can impose a priori con- 
straints on the field topology - e.g., that the field is purely 
potential, or purely toroidal, or a combination of both. An- 
other important advantage of this formalism is that both 
simple and c omplex magn etic topologies can easily be re- 
constructed (|Donatil l200iri . whereas the original method 
failed at reconstructing simple magnetic geometries (such 
as dipoles; Brown et al. 1991). 

To effect this approach, we describe the field as the sum 
of a potential and a toroidal component, each expressed as 
a spher ical-harmonic expan sion. In a formalism similar to 
that of IJardine et al.l (|l999h . the field components can be 
written as: 



Br{e,4>) ■ 
Be{e,<p) 

where 

Zi,m{0, (j)) 



^tt€,ml£,m(^,0) (2) 

i,m 

{Pi,mZe,m{0, 0) + 'Ji,mXe,m{0, 0)) (3) 

i,m 

{Pi,mXi,m{0, 0) - -fi,mZi,m{0, 0)) , (4) 



dPi,m{e) imcb 
i + 1 86 ^ 

Q,m Pi,m{0) . 

1 — 7 — — TT- I'm e 
sin 6* 



2^+1 {i-m)\ 
47T (^ + m)!' 



(5) 
(6) 

(7) 
(8) 



with i and m denoting the order and degree of the spherical- 
harmonic mode Yi^rn{0, (/)) {0 and being the colatitude and 
longitude at the surface of the star), and Pi,m{0) the asso- 
ciated Legendre polynomial. For a given set of the complex 
coefficients a^,m, f^i,m and ji^m (where ai^rn characterises 
the radial field component, Pi^m the azimuthal and merid- 
ional components of the potential field term, and 7^,rn the 
azimuthal and meridional components of the toroidal field 
term), one can produce the associated magnetic image at 
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the surface of the star, and thus derive the corresponding 
Stokes V dataset. We carry out the inverse problem, aimed 
at reconstructing a set of complex coefficients from an auto- 
mated, iterative fit to the observed circular-polarization LSD 
profiles. Principles of maximum-entropy image reconstruc- 
tion are applied to the set of complex coefficients, rather 
than on t he image pixels. Thi s is similar to what is pre- 
sented by Hussain et al. ("2001), except that we generalize 
the problem to fields that are non-potential and feature a 
significant toroidal component. Fitting a pure potential field 
to the data is equivalent to fitting and f3i^rn alone (set- 
ting all 7£,rn to zero); using all three sets of coefficients in the 
fitting procedure produces a more general magnetic topol- 
ogy, with a non-zero toroidal field. Trying both approaches 
is a straightforward way of investigating whether or not the 
magnetic field at the surface of r Sco is potential in nature. 

For a spectral resolution of 5 kms~^, a microturbulent 
velocity of 7.5 kms~^, and a projected equatorial velocity 
VeSin(i) of 6 kms~^, the number of resolved equatorial el- 
ements around the star is about 5, implying that we need 
about 10 equatorial elements to reproduce the observations 
with adequate surface sampling. Truncating the spherical- 
harmonic expansion of the magnetic-field components to 
terms with ^ < 10 is therefore sufficient in the case of r Sco, 
and should introduce negligible degradation in the spatial 
resolution of our reconstructed images. In practice, we used 
I < 12 and found that, as expected, no improvement in 
the quality of the fit to the data was obtained when adding 
higher-order terms. This corresponds to mapping a total of 
90 modes at the surface of the star, implying a total of 360 
image parameters in the case of a potential field, and 540 in 
the case of a potential- plus toroidal-field configuration. 

We first carried out a series of magnetic reconstructions 
for a wide range of values of the rotation period (without 
constraining the field to a specific type of configuration). 
The minimum is obtained at Prot = 41.02 zb 0.03 d (la 
uncertainty). This estimate of Prot is a refinement of that 
derived in Sec. 13.11 because of the more complete physical 
model; it is fully compatible with the adopted, more precise 
period derived independently from the lUE data (eqn. Q 
Sec.|S2J. 

The best-fit model of the Zeeman signatures is shown 
in Fig. [7| from which it is obvious that the greater part of 
the observed profile information is satisfactorily reproduced. 
Nonetheless, the minimum is as large as 1.5, indicating 
that discrepancies between the model and observations still 
exist. The origin of these (small) discrepancies is not yet 
clear, but may result from the simple isotropic local line- 
profile model we use to compute the synthetic Stokes V pro- 
files (Sec. 21- We note that the fit to the data is much worse 
if we force the magnetic topology to be very simple, e.g., 
similar to that found in most magnetic chemically peculiar 
stars to date. For instance, when truncating the spherical 
harmonics expansion to ^ < 1 (equivalent to fitting the data 
with a tilted magnetic dipole model), the minimum achiev- 
able is 15; with I <2 (roughly equivalent to adding up 
a magnetic quadrupole component to the model), the fit 
quality is still very rough (x^ — 8). The detected Stokes V 
profiles (and in particular the observed rotational modula- 
tion) definitely indicate that the magnetic-field topology of 
r Sco is much more complex than usual (by massive star 
standards) . 



5.2 Results 

The magnetic topology we reconstruct by assuming that the 
field includes both a potential field and a toroidal field com- 
ponent is shown in Fig.|Hl The corresponding x? is 1.5. If we 
instead assume that the field is purely potential, the optimal 
fit we obtain yields x? — 1-8, indicating a slightly poorer 
fit; while the inferred topology is broadly similar to that 
shown in Fig. |H] the contrast and information content is sig- 
nificantly higher than in the previous case. We can thus say 
that, at given image information content, the potential- plus 
toroidal-field configuration provides a significantly better fit 
than the pure potential-field configuration (with respective 
xi values of 1.5 and 1.8), indicating that the first option 
appears significantly more likely than the latter. 

The contributions of potential and toroidal terms to 
the azimuthal and meridional field components are shown 
separately in Fig. |^ (the toroidal term does not contribute 
to the radial field component). We therefore conclude that 
the surface magnetic topology of r Sco is mainly poten- 
tial, but apparently also includes a toroidal field compo- 
nent. We find that the potential field component includes 
about 70% of the overall reconstructed magnetic energy, and 
thus clearly dominates over the toroidal field component. In 
particular, we note that, compared to the potential compo- 
nent, this toroidal component is much weaker in r Sco than 
in most cool magnetic stars observed to date, where the 
toroi dal comp onent largely exceeded the poloidal one (e.g., 
Dona ti et al. fiooa). Moreover, we find that this toroidal 
component is not primarily axisymmetric (at least about 
the rotation axis), featuring for instance two sign switches 
(at phases 0.52 and 0.18) at equatorial latitudes (see Fig.|^ 
top right panel). 

The results further emphasize that the reconstructed 
magnetic field is far more complex than a simple dipole; for 
example, the polarity of the radial field component switches 
sign six times along the equator (instead of just twice as 
expected for a tilted dipole). A more quantitative way of 
considering the complexity of the field is to examine the 
relative strengths of the reconstructed spherical-harmonic 
coefficients a£,rn, /^£,m and 7£,m (see Fig. llOf) . Even though 
the modes corresponding to a tilted dipole are excited (^ = 1 
and m = — 1), a large number of other modes carry sig- 
nificant power (especially sectoral modes); the I — m — A 
mode is among the strongest, and its signature can be read- 
ily seen from the magnetic images of Fig. |H| Unsurprisingly, 
modes with I > Q contribute little to the reconstructed im- 
age. The reconstructed toroidal component only shows up 
in low-degree (though non- axisymmetric) modes. 

Since our r Sco spectropolarimetry spans several rota- 
tion cycles, and since several phase ranges were observed 
more than once (phases -0.00, 0.47 and 0.85; see Tabled, 
we can directly investigate whether the magnetic field of 
T Sco exhibits signs of variability on a timescale of about a 
year. We find that profiles obtained at very similar phases 
in different cycles (e.g., at cycles 2.02 and 10.01) agree to 
within the noise level. 

We can also estimate how much latitudinal shear the 
magnetic topology of r Sco experienced between 2004 Sep. 
and 2005 Sep., us ing; t he methods em ployed for cool stars by 
IPetit et al.1 j20o3) and lPonati et all |2003); we find no evi- 
dence for differential rotation at the surface of r Sco, with 
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Figure 7. Maximum-entropy fit (thick lines) to the observed Zeeman signatures of r Sco (thin hues). The rotational phase and cycle of 
each observation is written next to each profile. A 3a error bar is also plotted left to each profile. 



an upper limit of about 3 mradd"^ (i.e., at least 20 times 
smaller than for the Sun). This is in agreement with our 
finding that the rotational modulation of UV lines, presum- 
ably related to the magnetic topology (Sec. 17.3.1^ . is stable 
over timescales of decades (Sec. 13.21) . 

We emphasize that only moderate surface spatial reso- 
lution can be obtained for r Sco, as a result of its low rota- 
tion velocity; close bipolar groups, for example, could there- 
fore easily remain undetected if present on scales smaller 
than the resolution element. However, given the dense phase 
coverage obtained throughout the whole rotation cycle, the 
large-scale magnetic field of r Sco (up to orders with ~ 6) , 
as well as its non- variability on a time scale of about 1.5 yr, 
is very well constrained by our observations. 



6 ORIGIN OF THE FIELD 

Using the magnetic map we derived for r Sco, and thanks in 
particular to its unusual degree of complexity (by the stan- 
dards of early-B and 0-type stars), several questions can be 
addressed regarding the physics of massive stars. Our results 
both give us the opportunity to rediscuss the problem of the 
origin of magnetic fields in very hot stars, and also enable 
us to investigate the impact of complex fields on radiatively 
driven winds. 

Although the classical picture was that magnetic fields 
of hot stars were presumably fossil remnants from the 
formation stage, the situation has changed considerably, 
with regular reports from both observers and theoreti- 
cians that massive stars may be able to generate dy- 
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Figure 8. Maximum-entropy reconstructions of the magnetic 
topology of r Sco, assuming that the global field can be expressed 
as the sum of a potential field and a toroidal field. The three com- 
ponents of the field are displayed from top to bottom (fiux values 
labelled in G). The top image (radial field component) is de- 
scribed through the set of complex coefficients a^^rn (see Sec.|3. 
The star is shown in ffattened polar projection down to latitudes 



namo processes, either deep inside their convective cores 
irbonneau MacGregoJ l200ll : iMacDonald MullanI 
I2OO4I : iBrun et a h '2005*1 , within the greater part of their ra- 
diativ e envelope (Spruit 199 9, 2002; MacDonald & MullanI 
2004l:lMn11an X. MacDona1dlEoo3 iMaeder ^. MevnetllS 
Braithwai tell2006l) , or in a subsurface layer (|Tout Pringld 
1995; Lign ieres et al.|[l99^ . In each case, different processes 
are invoked to explain the generation of magnetic fields. How 
do these proposals stand up in the particular case of r Sco? 

6.1 Dynamo processes? 

Being intrinsically a very slow rotator (the s econd slowest 
rotator among massive stars after HD 191612; iDonati et al.l 
2006), T Sco does not seem to be an optimal candidate for 
triggering dynamo processes. Of course, one may argue that 
the rotation rate of the inner stellar regions may be far larger 
than that at the surface; evolutionary models tend, how- 
ever, to indicate that the radial gradient in rotation rate is 
only moderate in massive stars, particularly in those hosting 
magnetic fields (Maeder & MevnetJ i2003l . 12004 1200.^) . 

Wherever they operate (whether in the convective core, 
the radiative envelope, or a subsurface layer), dynamo pro- 
cesses are all expected to strengthen with rotation rate and 
to vanish when rotation is slow; they should therefore be 
relatively weak in a star like r Sco. I n the particu lar case of 
the Spruit-Tayler dynamo processes (ISDruitll999l. E002\ for 
example, field strengths only of order a few G are expected 
to appear at the surface of a star with a rotation rate as 
small as that of r Sco ( Mullan & MacDonald 2005), much 
lower than we have found. In the case o f the core-dynamo 
hypothesis, iMacDonald MullanI (12004) demonstrate that 
surface magnetic fields need to originate in core fields that 
largely exceed the equipartition value, which is again highly 
unlikely in a slowly rotating star. 

A second observation is that the magnetic features we 
reconstruct at the surface of r Sco are present at essentially 
all latitudes. This is not compatible with the predictions of 
the Spruit-Tayler dynamo model with buoyancy included 
to allow magnetic flux rise up to the surface, which should 
produc e magnetic regions concentra ted at intermediate lat- 
itudes (jMullan M acDonaldl l2005l) . Nor is it compatible 
with the core-dynamo theory, which predicts that flux tubes 
should mostly show up very close to the pole at the surface 
of the star (Ma cDonald k Mullanl|2004 l) . 

Our observations further indicate that the recon- 
structed magnetic field is mostly poloidal, and includes no 
more than a moderate surface toroidal component. This is 
not compatible with dynamos operating in a shear layer be- 
low the surface or within the convective zone, which are 
expe cted to produce s trong axisymmetric azimuthal fields 
re.g. . 'Brai thwaitell20 06^ that should likely show up at photo- 
spheric level, and should even dominate the global magnetic 
map in the case of a sub-surface shear layer (as they do 
in stars with very sh allow subsurface convective zones; e.g., 
iMarsden et al.lEoO^ . Moreover, we observe that the photo- 
sphere of T Sco experiences negligible latitudinal shear (at 
least 20 times smaller than that of the Sun). Again, this 
is most probably inconsistent with a subsurface shear-layer 
dynamo, which is expected to generate azimuthal and radial 
gradients of angular velocity, at least within the shear layer 
itself. 



12 J.-F. Donati et al 




Meridional magnetic field Meridional nnagnetic field 

0.50 0.50 




0.00 0.00 

Figure 9. Azimuthal and meridional components of the reconstructed potential (left column) and toroidal (right column) field structures. 
Adding both together yields the azimuthal and meridional field components shown in Fig. [s] The image on the left hand side is described 
through the set of complex coefficients Pi^rn while that on the right hand side is obtained through the coefficients Ji^rn (see Sec-I^. 



An additional argument against dynamo processes is 
that they should essentially succeed (if conceptually valid) 
at producing magnetic fields in most hot stars and not only 
in a small fraction of them. The fact that magnetic fields 
are detected in a star like r Sco, known for its peculiar spec- 
troscopic morphology (e.g., through its abnormally strong 
UV P-Cygni lines and its unusually hard X-ray emission), 
after having been detected in other peculiar hot stars (like 
0^ Ori C, HD 191612 and Cep), represents further evi- 
dence that magnetic fields (at least those of moderate to 
high intensity) are not a common feature of most hot stars, 
but rather a rare occurence. In this respect, the parallel 
between massive hot stars with intermediate-mass ones, for 
which only about 10% of them are magnetic, seems strength- 
ened by our new result. 



Finally, the magnetic configuration of r Sco exhibits no 



sign of variability evident on a timescale of a year or so.^ 
Archival lUE data provide a further, albeit indirect, indica- 
tion that the field topology is probably stable on timescales 
of decades. Although all dynamo models proposed for hot 
stars predict some kind of temporal variability, little infor- 
mation is available on the typical timescale on which the 
surface fields are expected to evolve. At least for models 
depending on differential rotation, we can expect the field 
configuration within the star to evolve on timescales shorter 



^ We note that two longitudinal-field estimates obtained for r Sco 
by Landstreet (1982), -61 ± 30 G and -16 ± 23 G and on JDs 
2,442,174.92 and 2,442,175.89 respectively, average to a longitu- 
dinal field of —33 ± 18 G at a phase 0.49 (cycle —169) accord- 
ing to eqn. The longitudinal field we measure at this phase, 
— 20± 3 G, is compatible with the estimate given bv iLandstreeti 
^982) to within the errors, and provides additional (though ad- 
mittedly fairly weak) evidence that the magnetic field of r Sco 
exhibits no intrinsic variability on timescales of decades. 
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Figure 10. Modulus (in G) of the spherical-harmonic complex coefficients ai^rn, f^e,m and 7£,rn (see Sec. |5j for the reconstructed 
magnetic-field topology of r Sco, as a function of mode degree i and order m. These sets of coefficients respectively correspond to the 
magnetic images shown in the top panel of Fig. |3 (radial field component), left panels of Fig. IHl f azimuthal and meridional component 
of potential field term) and right panels of Fig. lol fazimuthal and meridional component of toroidal field term). Only modes with ^ < 10 
are displayed here. 



than a year, as buoyant flux tubes typically need timescales 
of order a year to travel to the surface in a star like r Sco 
(provided the field in the flux tube is not too weak com- 
pared to the local equipartition value: iMacDonald MullanI 
l2004r ). again implying that the observed non- variability ar- 
gues against the proposed dynamo m odels. 

Recent numerical results by iBraithwaite Suruitl 

(|2004ri and IBraithwaite Nordhmdl J200(tl also indicate 
that magnetic fields in hot stars seem to reach a stable 
equilibrium (involving both potential and toroidal fields) 
but are not self- amplified by instability processes, unless the 
star features self-sustained differential rotation ( Braithwaite 
1200 6"). Similar numerical results are obtained by Brun and 
collaborators (Zahn, personal communication). Our result 
indicates that r Sco is likely not a differential rotator (as it 
would otherwise host strong axisymmetric toroidal fields like 
those seen on cool stars, Donati et al. 2003); it may there- 
fore be fairly natural that no evidence for dynamo action is 
detected on r Sco. 



6.2 Fossil fields? 

The next step is to compare our observations with predic- 
tions of the fossil-field theory. Since r Sco is rather young (a 
few Myr; Sec^J, the complexity of the field we detected is 
probably not a problem; while low-order terms (with longer 
decay times) are expected to dominate the fossil magnetic 
topologies of old stars, higher-order terms should still be 
present in stars as young as r Sco. Moreover, very little dif- 
ferential rotation and variability (on a ^yearly timescale) is 
expected to occur in stars hosting superequipartition fossil 
magnetic fields, in agreement with what we find. 

Both toroidal and poloidal fields of comparable 
strength are expected to be present within the star, at 
least to ensure dynara ical st ability of the fossil field on 
long timescales (IMossI EoOll: IBraithwaite Spruitl 120041: 
IBraithwaite Nordlundl l2006l) . The prediction is that the 
expected toroidal field should be roughly axisymmetric 
(with respect to the poloidal magnetic axis) and concentrate 



on the poloidal magnetic equator (|Braithwaite Spruitl 
12004 IBraithwaite Nordhmdl EoO^ . Although we indeed 
detect a small toroidal field at the surface of the star, its 
topology is not compatible with such predictions, which 
would require the toroidal field to coincide with the poloidal 
field equator (i.e. to show up mainly as a meridional field belt 
encircling the star, passing through both rotational poles 
and crossing the rotational equator at phases of about 0.3 
and 0.8, see Sec. |7|. Note however that theory expects the 
toroidal field to remain within the stellar interior; it may 
therefore be unsurprising not to detect it at photospheric 
level. We suggest that the toroidal field structure we detect 
at photospheric level rather results from the interaction of 
the stellar wind and the magnetic field (see Sec.0. 

Another attraction of the fossil-field theory is that there 
is no contradiction with the fact that the star is both mag- 
netic and slowly rotating. Actually, we note that magnetic 
hot stars are, in average, even more slowly rotating than 
non-magnetic stars. One can, of course, wonder whether this 
is a real property of massive magnetic stars, or simply an 
observational bias (magnetic fields being easier to detect by 
spectropolarimetric methods in narrow- lined, slowly rotat- 
ing stars); however, experiments clearly demonstrate that 
several rapidly rotating O stars (e.g., C Pup, C Ori) have 
surface magnetic fields with strengths not more than a few 
tens of G (Donati, in preparation). The youth of r Sco 
excludes the possibility that slow rotation is a result of 
angular- momentum loss through a magnetic wind during the 
main-sequence phase;^ thus the most probable option is that 
this situation is due to a process occuring during the for- 
mation stage. The idea proposed for magnetic, chemically- 
peculiar stars (also more slowly rotating in average than 



^ Following iDonati et all J2006I) . we evaluate the magnetic- 
braking timescale of r Sco to be of order of 5 Gyr. As this 
timescale is some 3 orders of magnitude larger than the age 
of r Sco, we can safely conclude that angular-momentum loss 
through the current magnetically confined wind is not responsi- 
ble for the slow rotation. 
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non-magnetic stars of similar spectral type) , invoking mag- 
netic coupling with a putative accretion disc (IStepienl 200(1 V 
is probably not applicable in the case of massive stars, which 
do not exist as stars during the formation stage and directly 
appear onto the main sequence. One possibility is that proto- 
stellar discs with intrinsically higher primordial magnetic 
fields are more successful at expelling angular momentum 
from the disc (e.g., through magnetic jets) than those with 
weak primordial fields, leading to magnetic hot stars rotat- 
ing more slowly than non-magnetic counterparts of similar 
mass. 



6.3 Conclusion 

Taking all these arguments into consideration, we find that 
the magnetic topology we have reconstructed for r Sco is 
more likely to be of fossil origin than to be generated by 
any of the various dynamo mechanisms proposed up to now 
in the literature. If this is confirmed, it would indicate that 
very hot magnetic stars probably represent a high-mass ex- 
tension of the classical Ap/Bp phenomenon; the reason these 
massive magnetic stars do not mark themselves as chem- 
ically peculiar is probably related to their strong winds, 
which prevent photospheric element stratification building 
up. This scenario would also argue i n favo u r of the proposi- 
tion of iFerrario Wickramasinghel (|2005l boOGf) , who sug- 
gested that massive magnetic stars are the progenitors of 
highly magnetic neutron stars. 



7 THE EXTENDED MAGNETOSPHERIC 
STRUCTURE 

A second topic of interest is the impact of the magnetic field 
on the radiatively driven wind of r Sco. In particular, r Sco 
gives us the opportunity of investigating the confining ef- 
fect of magnetic fields whose topology is more complex than 
those of other early B and O stars for which similar studies 
have been carried out (iDonati et al.l200lLEoo3 : lGagne eraP 
l2005b aV 

In the now-sta n dard picture, initially proposed by 
( ^97 ) and further investig at ed b y 

200211 ■ C^^^^^^L^OOl' 

Tbwnsend k Owocki (2005) and Gagne et al. (2005b a), the 
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magnetic field is assumed to be dipolar; the dense wind com- 
ing from each magnetic hemisphere is deflected by the field 
towards the magnetic equator, where it produces a strong 
shock, a very hot X-ray emitting post-shock region (reach- 
ing temperatures of 10^ K), and a cool, dense disk in the 
magnetic equator, where the plasma accumulates before be- 
ing ejected away from, or accreted back onto, the star (de- 
pending on the local radial velocity of the plasma when it 
reaches the disc, and on the effective gravity in the disc at 
this point). 

For a more complex magnetic topology, the picture 
is expected to differ significantly. The extended magnetic 
structure should show a correspondingly greater degree of 
complexity, involving distinct regions of closed loops alter- 
nating with regions of open field lines, rather than two open- 
field polar cones and one closed- field magnetic torus, as in 
the dipole field case. Wind flows should freely escape the star 
along open field regions and should produce shocks and very 



hot X-ray emitting plasma within each closed-field region, 
with cool, dense condensations forming at loop summits. 
The resulting magnetospheric structure should therefore be- 
gin to resemble that of the Sun, with hot coronal arcades 
confining cool, dense, prominence-like structures (the main 
difference being, of course, the heating mechanism itself). 



7.1 The magnetic-confinement parameter 

First, to evaluate to what distance from the star the wind 
of T Sco is magnetically confined, it is useful to consider 
the wind magnetic-confinement parameter 77, defined by 
lud-Dou~ fc Owocki (2002) to characterize the ratio between 
the magnetic- field energy density and the kinetic energy den- 
sity of the wind: 

r] = B^rJ/Mvoo, ^ (9) 

where B is the typical magnetic- field strength, M is the av- 
erage mass-loss rate, and Voo is the terminal wind velocity. 
Our study demonstrates that B is in average :^ 300 G over 
the surface of r Sco. However, there is a substantial disper- 
sion in published estimates of M, with r ecent determinations 
ranging 0.2-6 x 10"^ Moyr"^ (e.g., iMokiem et all 1200.4 
ReDolust et al. 2005); we adopt a value of 2 x 10~^ M© yr~^ 
as a reasonable average of modern observational determina- 
tions. Observation s also suggest a terminal velocity of about 
2,000 kms"^ (e.g.,|Abbot3E92i), but line-driven wind theo- 
ries predict a significantly higher value (up to 3,800 kms~^, 
Pauldrach 1987; note, however, that the stellar parameters 
used for these theoretical studies, and in particular the mass, 
temperature, luminosity and radius of r Sco, are all signifi- 
cantly overestimated compared to modern values). We adopt 
= 2, 000 kms"^ dMokiem et al.l l2005\ 
With these values, we find that the wind confinement 
parameter is about 40, and thus that the Alfven radius, 
above which all closed magnetic loops open under the wind 
ram pressure, is of order of 2 R^,. This is in good agreement 
with the findings of lCohen et al.l (|2003l) . based on Chandra 
observations (and in particular on Si Xlll and Mg XI line ra- 
tios), that the X-ray emitting plasma is located at average 
distances of about 1 R^ above the surface^^. 

Arguably, this agreement may be partly coincidental, 
given the large uncertainty on M; if it were 10 times smaller 
than the value we adopted, with all other parameters held 
fixed, ?7 would reach a value of 400. The magnetosphere 
would then be confined out to significantly larger distances 
(typically of order 4i?^); we should then detect hard X-ray 
emission from loops extending several R^, above the stellar 
surface, which is apparently not the case. If M were larger 
than 6 x 10~^ Moyr"^, then 77 would be smaller than 10, 
and no stable magnetic loops extending further than 0.3 i?^ 
above the surface would survive the wind pressure (e.g., 
lud-Doula Owockil boo'^) . again in contradiction with X- 
ray observations. We therefore conclude that, if this model 
is correct. X-ray observations constrain 77 to values ranging 



A recent re-analysis of the Chandra r Sco spectra using up- 
dated atomic data (e.g., iCagne et al.ll2005cJ) indicates that all 
line ratios from He-like ions (including both Si xill and Mg xi) 
are consistent with the X-ray emitting plasma being concentrated 
at a distance of about 1 Ri, above the photosphere (Cohen 2006, 
personal communication). 
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Figure 11. Closed magnetic- field lines of the extended magnetic configuration of r Sco, extrapolated from the photospheric map of 
Fig. 13 The star is shown at phases 0.25 (left) and 0.83 (right). Note the warp of the magnetic equator and the additional network closed 
loops around phase 0.65 (mostly visible on the right side of the right panel). 




Figure 12. Open field lines of the extrapolated extended mag- 
netic configuration. The star is shown at phase 0.83 only. 



typically within 20 to 100, and hence the mass-loss rate to 



values of 1-4x10 ^ M© yr 



7.2 



Magnetic- field extrapolation 

Using the field-extrapolation technique of Ijardine et alJ 
(|l999i) , with a spherical source function set to 2 (to mimic 
the mainly radial orientation of the field lines at distances 



larger than 2 i?^), we can draw inferences on the large-scale 
magnetospheric structure of r Sco, presented in Figs. ITTI 
(closed field lines) and 1121 fopen field lines). We find that 
the extended field structure is significantly more complex 
than a global dipole, even though it still features a torus 
of closed magnetic loops encircling the star (with an axis 
of symmetry roughly tilted at ^90deg to the rotation axis; 
Fig. ITT)) and two main cones of open field lines on opposite 
sides of the star (see Fig. I12f) . In particular, the magneto- 
spheric equator is significantly warped, and additional net- 
works of closed loops are present at low latitudes (e.g., one 
around phase 0.65, left panel of Fig. 1111 in conjunction with 
the equatorial region of positive radial field reconstructed at 
this phase, another one around phase 0.4). Most closed loops 
typically extend to a distance of up to 2 i?^, in reasonable 
agreement with constraints derived from Chandra data. 



Note that these small networks of closed loops (at 
phases 0.40 and 0.65) both coincide with local maxima of the 
reconstructed toroidal field component (see top right panel 
of Fig. El; from this apparent spatial correlation, we spec- 
ulate that the toroidal field component we detect may be 
produced through an interaction of the stellar wind and the 
magnetic field right at the photospheric level. We also note 
that these loops roughly coincide with the rotation phase at 
which the unpolarised line profile of r Sco is slightly nar- 
rower than average; again, this may reflect the particular 
wind configuration that results from this specific field con- 
figuration at the stellar surface. MHD simulations are of 
course needed to confirm whether this idea is realistic or 
not. 
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7.3 Observational implications 

7.3.1 Optical and UV diagnostics 

This model implies that excess absorption in UV lines should 
occur when the magnetospheric equatorial plane crosses 
the line of sight, a s in all simil ar hot magnetic st ars (e.g., 
iDonati et a'lll200lL Eo02. 20od: iNeiner et aLlEooil . i.e., at 
phases 0.3 and 0.8 - in excellent agreement with observa- 
tions (see Fig. 1^. In particular, the fact that the two UV 
absorption events we detect are largely similar in shape, and 
separated by 0.5 rotation cycles, provides independent con- 
firmation that either the angle i of the rotation axis to the 
line-of-sight (which we found to be c:^ 70°), or the global 
tilt of the large-scale magnetic structure to the rotation axis 
(for which we derived an estimate of ^ 90°), or both, is/are 
large. As mentioned previously, the relative phasing between 
our new spectropolarimetric data and the old archival lUE 
spectra is accurate enough (of order 1%; Sec. 13. 2|) to en- 
sure that this match is not a coincidence; it thus provides a 
strong argument in favour of the present model. 

Photometric measurements secured by Hipparcos in- 
dicate a constant flux level (to within 10 mmag), show- 
ing that the column density of the wind material trapped 
within the magnetospheric equator of r Sco is not high 
enough to produce detectable light variations through scat- 
tering, even when the disc is seen edge-on. This situation 
is similar to that found for f3 Cep and Ori C (for which 
no ph otometric variations are detected; iDonati et al.ll200lL 
120021) . but differs from HD 191612 and a Ori E (for which 
eclipses of the continuum radiation by the magnetospheric 
plasma are obse rved, at levels of 0.04 and 0.15 mag, re- 
spectivelv: Walborn et allEooi iTownsend Rr^ OwockilEooi 
IDonati et al]l200(tl . 

In our data, we detect no rotational modulation of 
the Ka flux from r Sco, mimicking what is observed in 
f3 Cep (which shows only long-term Ha variations), but 
very different to 6* Ori C, a Ori E an d HD 191612 (all 
of which exhibit strong Ha modulation: Jstahl et allll996l : 
iTownsend Owocki 2005VWalbor n et al.l20 03V More mod- 
elling is required to check whether these observations are 
compatible with the basic picture presented here; this is 
postponed for a future study. 

7.3.2 X-ray diagnostics 

The wind pressure at the base of the the postshock region 
is given by 

= e Mvoo/^TiR^^, (10) 

to first order, where e = 1/x^ — and x is the radial 

distance from the centre of the star at which the equilib- 
rium location of the shock front settles (in units of i?*). 
For a range of reasonable values of x ('^1.2-1.7), e remains 
roughly constant ('-^O. 10-0. 15), implying a wind pressure of 
about 20-30 gcm~^ s~^. This corresponds to a number den- 
sity of protons and electrons of about 10^° cm~^ for a post- 
shock temperature of order 10^ K within the loop. Again, 
this is in good agreemen t with the up per limits on the elec- 
tron density derived by ICohen et al.1 (|2003) from Chandra 
data. By assuming that the closed 'corona' (the equatorial 
magnetospheric torus and the small additional networks of 
closed loops) is filled with such a plasma, and using the 



simple coronal-structure model of Ijardine et al.l (|2002fl . we 
find that the resulting emission measure is of order of a 
few 10^^ cm~^, in reasonable agree ment with actual mea- 
surements from X-ray spectra (Woi dowski S chulz 20051). 
This result essentially indicates that, in the model we have 
devised, the observed X-ray emission can be mostly at- 
tributed to the magnetosphere (and that the model is there- 
fore broadly consistent with X-ray observations). 

If our speculation is correct, it implies that the X- 
ray emission of r Sco should be modulated on a timescale 
equal to the rotation period (i.e., 41 d), as a result of the 
magnetosphere being partially eclipsed by the stellar disc 
in appropriate viewing configurations. We estimate the ex- 
pected fractional modula tion to be about 40 %, i.e., compa- 
rable to that of e Ori C (iGagne et al.ll2005al lli3). but should 
feature two main eclipse episodes each rotation cycle, cen- 
tred on phases ^0.3 and 0.8 (concomitant with the UV line- 
absorption events witnessed in lUE spectra); this should be 
easily detectable given adequate temporal sampling. 

Checking such predictions should provide a strong test 
of our model. If the observed X-ray rotational modulation 
is much weaker than expected, it could imply that most of 
the magnetospheric emission is produced in small-scale loops 
evenly spread over the stellar surface; this would argue for 
additional, high-order components of the magnetic topology 
that we are not able to detect in this study (as a result of 
the Umited spatial resolution provided by Doppler imaging 
for stars rotating as slowly as r Sco). 

At this stage, a more accurate model is obviously nec- 
essary to confirm the conclusions of the present paper, and 
to develop in more detail how well the X-ray spectrum of 
r Sco can be reproduced once the magnetospheric structure, 
including the wind-induced expansion of magnetic loops, is 
consistently taken into account. Such a sophisticated model 
should also aim to reproducing the observed modulation of 
UV spectra and the upper limits on the pho tometric and Ha 
flux variability (e.g., as done for a Ori E bv ITownsend et al] 
l200.^ . 



8 CONCLUSION 

We have reported the detection of a magnetic field on the 
massive BO. 2 V star r Sco, using data obtained mostly 
with ESPaDOnS, the new high-resolution stellar spectro- 
polarimeter recently installed at CFHT. From the Zeeman 
signatures and their temporal variability, we were able to 
identify the rotation period of r Sco and to reconstruct the 
large-scale topology of its photospheric field. Archival lUE 
spectra confirm that the rotational modulation is stable on 
timescales of decades, with Prot — 41.03 d. 

We find that the surface magnetic topology is unusually 
complex (judged by the small sample of massive- star results) 
and is mostly potential. It also includes a moderate toroidal 
component; in particular, the strength of this toroidal com- 
ponent (relative to that of the poloidal component) is much 
lower than that found in partly-convective cool stars hosting 
dynamo- generated magnetic fields. No temporal variability 
of the magnetic structure is detected over the 1.5-yr period 
of our observations; we thus conclude that any surface dif- 
ferential rotation of r Sco is at least 20 times weaker than 
that of the Sun. 
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We determine that the large-scale magnetospheric 
structure of r Sco is significantly more complex than a global 
dipole; it features in particular a significantly warped torus 
of closed magnetic loops encircling the star, tilted at about 
90° to the rotation axis, as well as additional (smaller) net- 
works of closed field lines. The extended magnetic topology 
we derive from extrapolations of the photospheric magnetic 
maps is apparently compatible with the published X-ray lu- 
minosity and spectral characteristics of r Sco. Our model is 
also compatible with the observed modulation of UV spec- 
tral lines. We predict that r Sco should exhibit a clear ro- 
tational modulation of its X-ray emission. 

From these results, we conclude that its magnetic field is 
most probably a fossil remnant from the formation stage. We 
cannot yet completely rule out the possiblity that the field 
is produced through one of the recently-elaborated dynamo 
processes that may operate in the radiative zones of hot 
stars, but our findings already indicate that this option is 
rather unlikely. 
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